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I. Summary 

At the Bartol Coude Observatory we have continued our 
observavicns of umbral cores, both by multicolor photometry 
and by narrow -band photometry in the vicinity of the sodium 
D lines. These observations have indicated that no single 
umbral model can account for all cf the observed umbral cores. 

To speak of "the" model of an umbra is therefore net meaningful. 
In a series of observations cf one spot at a variety of posi- 
tions on the disk, it has been found that the sunspot exhibits 
limb brightening. Furthermore, the weakening of sodium ab- 
sorption near the limb (relative to that at disk center) is 
much more pronounced than any model currently predicts. Both 
these results imply that the upper atmosphere of this spot is 
hotter than any current model predicts, and may in fact have 
a temperature inversion at the level where the D lines are 
formed. Of the models currently available we believe that 
Kneer's (1972) model comes closest to agreeing with the obser- 
vations, although even that model cannot explain all of the 
available observations. 

Theoretical work on Alfven wave propagation in the atmos- 
phere of a spot has been followed numerically to first order 
in the wave amplitude. In regions where the wavelength exceeds 
the density scale height, the wave flux is found to decrease with 
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altitrde in accordance with the WKB solution. But in the upper 
photosphere and chromosphere, where the density scale height 
is small compared with the wavelength, the first crder solution 
indicates a flux diverging at great height. Attempts to 
incorporate second-order terms in the wave amplitude are con- 
tinuing. 

Our observations of excess flux around sunspcts suggest 
that a significant fraction of the missing flux in sunspots 
may appear as a slight brightening in a large area surrounding 
spots. A time -dependent model of turbulent heat transfer 
indicates that turbulent diffusion by supergranules can account 
for the observed excesses in some cases. Sound waves also cai 
transport energy laterally a way from spots. 

The Alfven-wave model developed for sunspots has been 
extended to the case of starspots on flare stars and cn the 
primary component of contact binaries. The model has been 
successful in accounting for a wide variety of observations 
in these objects, and this gives us further confidence that 
the models developed for sunspots are indeed realistic. 

Effects of magnetic fields on solar neutrinos are shown to 
be negligible. 
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II. Observational Work with the Bar to I Coude Telescope 

a. Narrow Band Observations of Jxbral Cores 

Observations of the intensity of sunspots have been con- 
tinued on those clear days when suitably large spots were 
present on the disk. Solar activity has unfortunately declined 
during the past year, and we have observed only about a dozen 
spots this year, costly of small umbra 1 area. 

IiS observations which we have concentrated cn are an 
investigation of the maximum intensity between the sodium D- 
lines. In figure 1, we shew an idealized porticn of the spec- 
trum of a spot in the vicinity of the D- lines . We measure 
intensity, both in the umbra 1 core and at the center of the 
solar disk at three wavelengths, 5874.125 A, 5893,4 A, and 
5931.75 A, using wavelength resolution of 0.2 A. The first 
and third wavelengths are at the center of continuum windows 

tf 

in the umbral spectrum, as listed by Wohl (1970). Interpola- 
tion in wavelength between these two windows allows a deter- 
mination of the continuum intensity 1^ at 5893.4 A, which is 

the wavelength where the umbral spectrum has maximum intensity 

«« 

between the D- lines according to the Gcttingen umbral atlas. 

1^ is obtained beth in the spot and at the solar center. 

The direct measurements of the intensity at 5893.4 A yields 
1 B again both in the spot and at the center of the disk, In 
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the disk spectrum, where the D- lines are narrow enough that 
the wings dc not overlap appreciably, (sun) » I A (sun). 
Me therefore obtain 


‘spot 


Ig (spc-) m Ib (spot ) / 1 A vapor ) 


I A (spot) (sun) 


/< 


(sun) 


We use this, in combination with the intensity ratio 

l A (spot) 

R spot = i (sun) 


to define a two-dimensional classification system for umbra 1 
spectra. As with all umbral measurements, both r and R must 
be corrected for scattered light. We have assumed, as before, 
that by sitting in the umbral cere, where the spot intensity 
consistently falls to the same low level, we have essentially 
thereby Isolated the intensity of the umbral core, unaffected 
by blurring, and affected only by large angle scattering. 

The latter has been found to be a remarkably constant function 
for our instrument. The procedure cf correcting for scattered 
light is therefore simply to measure the aureole intensity at 
a distance of one arc minute off the limb of the sun, and 
multiply this by a factor Z. The value of Z becomes larger 
the closer the spot is to the disk center, and can be as large 
as 7-8. In presenting the results for spots near the disk 
center (Figure 2) we have indicated a range of values adopted 
for Z typically, Z * 4 to Z * 8. Near the limb, we choose 
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Z * 1-3 typically. The advantage of our two-dimensional 
classification system for unbral spectra is then immediately 
apparent, since scattered light corrections enter in a cor' 
related fashion in the two dimensions. Her.ce instead of 
uncertainties due to scattered light making their appearance 
as error boxes, in this diagram, the uncertainties appear as 
lines . The length of a line in the Figure is an indication of 
the scattered light in the instrument. On and before July 7, 
1973, the scattered light was almost 1% at one arc minute 
from the limb, causing the lines of July 7 and February 28, 

1973 to be very long. On and after August 30, 1973, follow- 
ing a cleaning of the sidercstat and 8-inch diaphragm, the 
scattered light fell to 0.4-0. 5%, and this is reflected in 
the shorter lines fcr the later observations. 

Theoretically, each sunspot model yields just one point 
in this diagram, at least for spots at disk center, p. * 1. 

The values of p. for different spots is indicated in the figure. 
Five models are shown, due to Kneer (1972), Stellmacher and 
Wiehr (1972), Yun (1971), and two revised models of Yun's 
«, reported in the annual report for the period September 72 - 
September 73, Section III. A), denoted YI and YII in the figure. 
All theoretical points move upwards and to the right as a spot 
moves towards the limb (but not enough to explain the results 
in Figure 2 , see below) . 
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Of special interest in Figure 2 is the series of observa- 
tions for a spot from July 9 to July 17, 1974, during which 
time the spot longitude varied from 64° East to 4G° West. The 
intensity of the continuum decreased near the center of the 
disk, indicating limb brightening from about 107o at disk center 
to 137» at 64° longitude. Note that all cur intensities are 
quoted relative to the intensity of the center cf the solar 
disk. If we were to give intensities relative to the solar 
intensity at equal distance from the disk center, the limb 
brightening would become even more pronounced. Thus 
I*C0-0)/l e {e-0) * 0.10, and I*(0*64°)/l®(0-64°) = 0.20. The 
good agreement between observations on July 11 and July 17 
when the spot was equidistant from the central meridian is 
an indication that the spot structure remained essentially 
unchanged in this six “day period. This spot therefore seems 
to be a good candidate for testing limb -darkening predictions. 
Since the observations at disk center (July 12) coincide well 
with Kneer' s model, it is encouraging to see that Kneer's 
model predicts limb brightening, although not as much as 307 o 
from n*l to m=0.6^ s (cf. Fig. 4 of Zwaan (1974) which is drawn 
for 7v= >000 £) . The model of Yun (1971) also predicts limb 
brightening at 5890 X, by a factor of 30%, although the contin- 
uum is too bright to be in agreement with present observations. 
On the contrary, the models of Stellmacher-Wiehr (1972) and of 
Zwaan (1974) predict limb darkening . The conclusion is that 
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for the umbral core in the spot which crossed the solar disk 
in 1974, July 9-17, the model which best fits the observa- 
tions is probably that of Kneer . 

Apart from the limb brightening in the continuum, the 
most obvious feature of Figure 2 is the fact that the line 
absorption is very much reduced near the limb. This is obvious 
not only from the observations on one spot in July but in all 
spots with |i-cos 0*0.6 or less. In all such cases, the line 
absorption is so weak that r spot *ises to 0,9, compared with 
0.5 - 0.7 near disk center. This is a sensitive test of 
umbral models, because it depends on two widely separated 
levels of the umbral atmosphere (t^O. 1 - 0.01 and t * 1-2). 
Stellmacher end Wiehr (1972) have computed the center- to- limb 
variation of r g p ot f° r 4 models (cf. their figure 2). From 
|i*l to |i*0.5, the value of r spot increases from 0.46-0.55, 
0.22-0.28, 0.58-0.73, and 0.77 to 0.80. Clearly none of 
these four models is adequate to explain the much larger 
dependence of r g p Qt on ij. which we have observed. All of the 
models give lines which are too deep near the limb, indicating 
that the temperatures in the upper layers of ail these models 
are too low. This conclusion is of course strengthened by 
the observations of limb brightening in the spot of July 1974. 
In fact these latter observations require the presence cf a 
temperature inversion, the amplitude of which cannot, however. 
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bc estimated without extensive computer modelling of the rad- 
iative transfer. In view of the sparseness of the currently 
available data, we believe that to embark on an extensive 
series of trial and error radiative solutions at this time 
would not be fruitful. Zwa&n (1968) points out the inevita- 
bility of a temperature inversion due to lateral influx of 
radiation. The inversion, he expects, should set in at ta.0.01 
in large spots, and perhaps as deep as t%0.1 in small umbrae. 
He presented some tentative evidence based on molecular lines 
to prove the existence of the inversion, although a synthetic 
spectrum analysis of the MgH molecular bands using the tempers - 
ture-inverted model has been shown by Branch (1970) to lead 
to disagreement with observations. Here we believe that we 
have obtained for the first time some evidence (at this time 
qualitative) based on atomic lines for the existence of the 
inversion in the atmospheric layers where the sodium lines 
are formed. 

This however does not change our conclusion that no 
single model fits all of the spot observations, and so every 
spot appears to be different, even in the deeper layers where 
the continuum is formed. 

In the continuum at 5893.4 A, the intensities of the 
umbral cores (R S p 0t ) extend over a wide range even near the 
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center of the disk, R spot ■ 6.5 - 24%. The lower limit which 
we have observed coincides well with the darkest spot observed 
by Ekmann and Maltby (1974) at > * 5790 A (I - 6.5 + 1.0%) while 
the upper limit which we have observed is much higher than other 
observers have indicated (cf. Figure 1 of Zwaan, 1974). 

It is apparent that (i) no two spots have the same physi- 
cal characteristics in both the atmospheric region where the 
wings of the sodium D lines are formed and in the region where 
the continuum is formed; (ii) none of the theoretical model 
points lies on any of the observed lines, suggesting that none 
of these models gives a satisfactory explanation of both line- 
forming and continuum - forming layers of a spec. In fact, 
the observations indicate that it is hardly meaningful to 
talk of ’'the 1 ' model of a sunspot umbra. Every umbra appears 
to be unique, at least within the ’-ange of umbral sizes in- 
cluded in the diagram, Ay * 8 to 60 millionths of a solar 
hemisphere. It is therefore difficult for us to agree with 

ft 

the suggestion of Rossbach and Schroter (1969) that umbrae 
of all sizes have essentially identical physical conditions. 
Although this might be the case in the deeper layers of the 
umbra (Henoux, 1969), there must arise differences in the 
upper layers, close to and above t *» 1 , where on the one hand 
photosphevic radiation begins to illuminate the dark gas in 
the: umbra, and on the other hand Lhe magnetic field geometry 
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i8 modifying itself in order to adjust to hydrostatic equilib- 
rium with the much lower gas pressures in the upper layers. 
Thus the assumption of a simple uniform field is acceptable 
at great depths, where it is permissible co write 


Pspct + B l' 8 ' * P sun 


( 1 ) 


(p = gas pressure* 5 - field strength). However, equation (1) 
is certainly not applicable within 200-300 km of the photosphere, 
where the pressure of a spot field of 2~3 kilogauss actually 
exceeds p sun » so that the above relation becomes meaningless. 

The field geometry presumably changes to a force- free config- 
uration by the time it reaches the chromosphere, and this 
change presumably begins at depths of 200-300 km below the 
photosphere, where curvature terms must be added to the left- 
hand side of equation (1). Since each spot finds itself in 
a different active region, each spot must adjust its deep 
field to different boundary conditions in the force- free 
region. Thus it is expected that the variation of magnetic 
field with depth is different for each spot, so that each 
3 pot has its own pressure variation, and accompanying T * T(t). 
The variations from spot to spot in Figure 2 may therefore be 
correlated either with the strength and geometry of magnetic 
fields in the active regions surrounding each spot, or with 
the diameter of the umbra, and its effect on lateral radiative 



influx. In the absence of detailed magnetographic coverage 
of each active region, we cannot at the moment test this hypo- 
thesis, but the need for magnetographic information in the 
vicinity of spots will also become apparent in discussing the 
excess flux around spots (see below). 

b. Broad-Band Observations of hinbral cores 

Two sjots were observed in broad band. One of these, 
observed on two days in July 1974, gave results as shown in 
Figure 3. The observations were corrects ’ for scattered light 
by applying a correction factor Z * 6 to the observed aureole 
intensity at one arc -minute from the limb. In Figure 3, the 
observations are compared with the five theoretical models 
referred to in Figure 2. On the basis of these broad-band 
observations, it appears that none of the five models is really 
suitable for this spot. From the results in Figure 3, it is 
difficult to avoid the conclusion that there are great differ- 
ences in intensity from one spot to another. The large spot 
(Ay « 80) observed on October 30, 1972, and plotted also in 
Figure 3 is clearly darker, by a factor of at least two, than 
the small spot (A u » 35). 

The efficacy of broad-band observations of spots is 
clearly minimal in view of corrections due to differential 
line blanketing which must be applied to the observations. 
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Our observations suggest that the variation of intensity cf 
umbra 1 cores with wavelength is net even monotonic . If this 
is due to errors in the differential blanketing corrections, 
then clearly breed-band observations cf spots cannot be used 
to set constraints on umbra 1 models, 

c. Magnetometer 

On July 17, 1974, a spot at u*0.76 was Investigated with 
the magnetometer. The observations were noisy, and they shewed 
circular polarization rising to +3.6X at L\ » 0.26 A. from the 
center of the line, however, the sign of polarization 
did not change going through the center of the line, but re- 
mained positive at AX both positive and negative. This may 
indicate that the field lines of the spot (viewed at an angle 
of 40° to the local vertical) were more or less transverse tc 
the line cf sight so that the. signal detected was actually a 
component cf elliptical polarization. 

d. Search for the Missing Flux of Sunspots . 

Following a suggestion of R.E. Danielson, we have under- 
taken an observational program which is aimed at searching the 
vicinity of sunspots for the possible presence of missing flux. 
It is possible that the flux which is missing from a spot is 
redistributed over a wide area surrounding the spot, sc that 
the intensity of the photosphere would be increased above 
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normal in the vicinity of the spot. Up to the present time, 
this excess of flux has never been detected around spots, 
partly because of limits of sensitivity of photographic 
techniques. In the worst cases , photographs might fail to 
detect excesses smaller than about 4% of the normal disk 
intensity (Bray and Lo^gfchead, 1964), sc that even if the 
missing flux were confined within an area of only 25 times 
the umbra 1 area (assuming for simplicity that the umbra is 
perfectly dark, and nc flux is missing from the penumbra), 
the excess flux would have been missed. Since the penumbra 1 
area is typically 4-5 times the umbra 1 area, the excess flux 
would have to be confined cc an area centered on the umbra 
with a diameter of nc more than 2.2 - 2.5 times the penumbral 
diameter, in order to be detected. (Ifae "bright -ring" obser- 
ved by Waldmeier (see Bray and Lcughhead , 1964; is a differ- 
ent effect, closely confined to the penumbral- photosphere 
boundary, and in all probability simply due to overshooting 
penumbral bright filaments (Danielson, 1973)). 

The fact that no observations have yet detected an excess 
flux around spots suggests that if this is indeed the correct 
explanation of the missing flux, then the excess must be spread 
out ever a much larger area of the solar surface, producing 
thereby an excess too small to be detected by photographic 
sensitivity. (However, very recently. Miller (1974) has 




succeeded in detecting the excess photographically in some 
cases. See below.) 

The sensitivity of the Bartel coude telescope is, however, 
essentially limited only by the electronics which prevent 
photon counting rates from exceeding 1G^ counts /second. Thus 
sensitivity of 0,1% is easily achieved, and if systematic 
effects can be excluded, the Bartel instrument is capable of 
detecting the excess flux if it were to be spread ever an area 
with diameter as large as 15 times the penumbra 1 diameter. 

We have searched for the flux surrounding spots on October 
26, 1973, (Ay * 15), November 15, 1973 (Ay =4), and January '3, 
1974 (Ay = 35). On the latter two days, sky conditions were 
not good enough over a lcng period to obtain spectral obser- 
vations of the spets concerned. In particular, on January 13, 
there were intermittent clouds, so that the results are hardly 
meaningful. The spot on November 15 was very small, so we 
will confine our remarks here to the results of October 26. 

The spot concerned was large, A ■ 370, but there was mere 
than one umbra, the largest having as umbra! area of only 
Ay * 15. The spot's position was 17°E 15°S (p. * 0.955). 

The telescope was stepped in declination starting at the 
center of the largest umbra, in steps of 2" per second, and 
the photon counts were recorded. Results are in Figure 4, 
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where each dot Indicates the counts in one second of time, 
starting in the umbra at radial distance r » 0. The scan 
continues to a radial distance of about 5 times the penombral 
radius, so the equivalent area is 25 times the penumbral area 
or 550 times the area of the largest umbra , 

The statistical uncertainty in the counting rates is 
indicated, and is obviously much smaller than the observed 
variations from second to second. The latter are probably due 
to real fluctuations in brightness in the photosphere, since 
cur entrance aperture, ll*6 corresponds very closely to the 
mean diameter cf granules (1100 km) on the solar disk. Granu- 
lar contrast is known to be as large as + 5% (Bahng and 
Schwarzschild, 1961), and this would certainly suffice tc 
explain the variations. 

More significant, however, is the gradual decrease in 
the counts with increasing distance from the spot. Ic is 
necessary first to remove limb darkening effects, for the scan 
extends along 10% of the solar diameter. The limb darkening 
has been calculated approximately along the scanning path, 
and is shown in Figure 4. It actually slopes upward towards 
the right, and so the gradual decrease of intensity away from 
the spot is enhanced. The dashed line through the points has 
an area in excess of the normal solar flux (represented by the 
curve labelled "limb darkening") which is sufficient to ccnpen- 
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sate fcr the area missing in Che spot, i.e. the area enclosed 
by the dash-dot line within the penumbra. 

Because the jog mode of guiding the telescope is most 
relia le In the north-south direction, we restricted cur scans 
t .■ flat direction. Thus, we have only a one-dimensional scan 
of what must surely be a two-dimensional feature. Nevertheless, 
it appears that we may indeed have succeeded in detecting the 
nisei ig flux around this spot, making its appearance in the 
form >f an excess of flux spread over a wide area surrounding 
the spot. We return to a theoretical discussion of the effect 
below. 

An interesting confirmation of our results 
has been provided by Miller (1974), who 

use \ photographic means to detect excess brightness around spots. 
The excess brightness observed by Miller may sometimes be as 
high ai 7 % of the normal disk intensity, and he believes that 
this excess may account for as much as 50% (if not more) of 
the mis. "ling flux. 

III. I ^oretical Studies at the Bartol Observatory 
a. Theoretical Zeeman Profiles of the Sodium T >2 Line 

while Dr. Kong Sik Yun was at Bartol, he had partially 
c upleted a computer program aimed at integrating numerically 
the equations cf transfer in the Stokes parameters I, Q, aid V. 
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This program has now been completed, and debugged sc that now 
the Zeeman profile can be computed for an arbitrary line in an 
arbitrary model atmosphere. In particular, the sunspot model 
of Kneer (1972) has been used to fit a theoretical profile to 
the observations of the D 2 line in a spot cn August 30, 1973. 
(The observations were reported in last year's annual report.) 
The resulting curve fits the observed points somewhat better 
than the poirts are fitted by theoretical results computed 
using Yun's spot model but the error bars on the observations 
are so large that there is little reason to prefer one model 
over the other on the basis of these observations. 

It has now become clear why Yun's initial estimates of 
expected circular polarization in the wings of the D lines 
were so small. Ke used Onno's (1956) theory for the Zeeman 
effect, which requires a knowledge cf the ratio of line to 
continuum absorption coefficients q at each point in the line 
profile. The procedure is to evaluate this ratio at the line 
center, where q = r> 0 , by fitting the observed line depth at 
line center with a theoretical line profile in the absence 
of any magnetic field. Then q can be evaluated at each point 
in the line profile by assuming that the line absorption 
coefficient has a Doppler profile, appropriately shifted in 
each polarization by the Zeeman splitting. However, in the 
program written by Yun, q was evaluated independently at each 
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point in the profile, and in particular was evaluated with 
physical conditions pertaining to one particular optical 
depth in the atmosphere, r » 0,001. Slice the effective 
level of line formation varies across the line, in the sense 
that in the far wings, the line is formed deeper than t 91 0.001, 
and in the core the line is formed at shallower layers than 
t * 0.001, the result is that the polarization is overestimated 
in the wings and underestimated in the core of the line. 

The program for numerical integration of the equations of 
transfer in polarized light has been modified to include more 
exact boundary conditions at the lower boundary of the atmos- 
phere. Rather than simply setting I = B (where B is the 
Planck function), and Q ■ V «= 0 there, we now use the more 
exact conditions suggested by Beckers (1969). However, the 
line profiles do not appear to be significantly affected by 
this change. 

A significant check on the numerical integrations is that 
the intensity I at the top of the atmosphere should approximate 
closely that evaluated by an eight- point formula for the inten- 
sity integral given by Zwaan (1965). If the numerical integra- 
tion has been terminated at an optical depth which is too large, 
then the numerically integrated line will not appear as deep 
as that predicted by Zwaan' s formula. We have found that at 
a distance of 0.6 A from line center, integration up to r * 0.001 
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en8ures accuracy within 0.67*. At larger distances from line 
center, it is net necessary tc carry the integration upwards 
as far into the upper atmosphere, while closer to the line 
center, in order to achieve accuracies of better than 17*, it 
is necessary to integrate higher up than t * 0.001. 

b. Zeeman Profiles of Other Lines 

In a preliminary attempt to extend the program to Zeeman 
profiles observed in magnetic stars, we have considered the 
profile of the line due to neutral chromium at 4254 A, a line 
which in the past has been observed in spots while liRSIES was 
located at the Flower and Cook Observatory in Paoli. This 
line is particularly interesting because it is the line which 
was used by Borra et al (1973) to search for weak magnetic 
fields in stars, using a Zeeman pclarimeter essentially iden- 
tical tc the one we use at Bartol for sunspots. Borra et al^ 
give a line profile observed in the magnetic star # Corcnae 
Borealis. We have computed Zeeman intensities of the Cr 1 
line using three model atmospheres for this star. Toe models 
were provided (by private communications) by S. Adel^an, J. 

Hardorp, and R. Friedman, with effective temperatures of 
4 n 

respectively 10 , 8750, and 8600 K, The peak in the observed 
polarization at A* * 0.12 A is reproduced, with a degree of 
polarization of 24%. However, these numerical integrations 
of intensity are only the first step in computing stellar 
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profiles, for these require further integration over the surface 
of the star. Such an integration would be extremely time- 
consuming, and hardly justifiable in view of our lack of 
detailed information about the distribution of magnetic fields, 
limb darkening, and microturbulent velocities over the surface 
of the star. 


c, Turbulent Heat Conduction and the Missing Flux in Sunspots 
Again following a suggestion of R.E. Danielson, we have 
investigated the horizontal transfer of heat away from a spot 
by turbulent convection. Wa idealize the spot as a cylinder 
j of radius a (» 10^* km) from which e: jrgy is to be transported 

j across the curved surface. The heat flux emerging from the 

•i 

! cylindrical walls of the spot is assumed to be equal to the 

29 

missing flux integrated over the total spot area, Q(«»I0 

3 

ergs/sec). The coefficient of thermal diffusivity outside 
the spot is set equal to a * 0.1 C g ax, where C g is the local 
sound speed, and * * ax is the "mean free path" of the turbu- 
lent conducting elements, expressed as a fraction x of the 
spot radius. The problem then is the following: given a 

region bounded on the inside by a cylinder of radius a, across 
which heat flows outward at a rate Q, what value must x (*\/a) 
have in order that the heat distribution outside the cylinder 
is compatible with our observations of the excess flux in tie 
neighborhood of sunspots? 
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To answer this, we need the solution to the cylindrical 
heat transfer problem, as tabulated numerically by Ingersoll 
et al (1950). Outside a cylinder transferring heat at a con- 
stant rate of q ergs per second per unit length of cylinder, 
the excess temperatures at radial distance r at time t is 
given by 

AT * f G(z, p) 
k 

where p » r/a, k * thermal conductivity, z * «t/a , a * k/pC p 
is the thermal diffusivity, p = gas density, Cp * specific 
heat at constant pressure. The function G(z,p) is an integral 

i. 

over zero and first order Bessel functions and Neumann ' s 
Bessel functions of the second kind. G is tabulated for four 
radial distances by Ingersoll et al, p * 1,2,5,10. (Note: 
p * 1 is the outside boundary of the cylinder) . The tabulated 
values are shown in Figure 5. 

For the purposes of computation, a model of the solar 
convection zone due to Baker and Temesvary (1966) was used. 

This gives values of C 8 , p, C p , and temperature as a function 
of depth below the top of the convection zone, i.e. as a function 
of depth below the undisturbed photosphere, essentially (since 
convection in the sun sets in close to optical depth t * 2/3; 
cf. Mu 1 lan, 1971). A value must be assumed for the depth of 
the spot, since AT above is expressed in heat flow q per unit 
length. We have pointed out reasons for taking the depth of a 
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spot equal to 10^ km (Mullan, 1972; 1973a) so this is what we 
assume here. The missing flux is taken to be 

Q - *a 2 F g (AyCl-Iy) + V 1 " 1 ?}) 
where F s = 6.33 x 10 10 ergs/cm 2 /sec is the solar flux; A. u , A p 
are the fractions of the total spot area occupied by umbra and 
penumbra respectively; and I u , I p are the intensities of umbra 

9 

and penumbra in units of the disk intensity. With a = 10 cm. 

Ay * 0.25, I u * 0.1, Ip * 0.7, we have Q * 8.95 x 10 2 ® ergs/sec. 

The problem is time dependent, so the bright ring around 
a spot diffuses outwards with time, and becomes weaker in inten- 
sity at any given point. The weakest the intensity can be 
occurs when the spot is as old as possible. Let us take a 
time t * 0.1 Ajjjgjj, which is the mean lifetime, in days, of a 
spot or spot group having a maximum area of A Tnay millionths 
of a solar hemisphere, (Bray and Loughhead, 1964). Then we 
can ask: what is the excess flux expected at the outer edge 
of the penumbra (p = 1) for a spot of area A max * 100 at age 
t * 10 days, given that x may vary over a range of values? 

We have set x * 0.003, 0.06, 0.03, 0.1 and 0.3 and calculated 
excess temperatures AT as a function of depth. AT is a rapidly 
decreasing function of depth, so that we need only consider 
the top-most L ; rs of the sun. As an extreme case, let us 
consider the level of the Wilson depression, which is about 


rn;: n 
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1000 km below the photosphere (Yun, 1970) . At this depth, the 
excess temperatures are found to be the following: 

TABLE 1 


X 

0.003 

0.01 

0.03 

0.1 

0.3 

AT(°K) 

4342 

2082 

1004 

425 

184 

AF/F(%) 

114 

55 

26 

10 

5 


The excess flux at the edge of the penumbra, AF, is computed in 

terms of the normal solar flux approximately by AF/F * 4AT/T 

where T is the undisturbed temperature in the sun at a depth of 

1000 km, T * 15242°K. It is apparent that unless x > 0.3, the 

excess flux even at this late time in the history of the spot, 

would have been easily detected even by photographic photometry. 

Is it possible to have x > 0.3? Presumably what is doing the 

turbulent convection around a spot is the convection zone, with 

diffusivity having such a value that supergranules can survive 

the erosion of turbulent conditions for only about 20 hours. 

With length scales of (2-4) x 10^ km, such lifetimes eorres- 

14 2 

pond to diffusivities of a « (0.6-2. 2) x 10 cm /sec. At 

A 

t - Kr km depth, C g * 12 km/sec, and hence such <*'s correspond 
to x * 0.5-2. Thus it appears that it might be marginally 
possible to explain the excess fluxes at the very end of spot 
lifetimes by this method. However, at earlier phases of a 
spot's life, say after one or two days, AF/F is expected to 
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exceed 5% even if x Is as large as 0.5-2, simply because the 
excess heat is crowded into a smaller area surrounding the 
spot. Further investigations are called for before a definite 
answer can be given, but, these preliminary results suggest 
that turbulent heat conduction away from spots is not really 
capable of explaining the missing f*.ux in spots. 

d. Sound Waves and the Missing Flux in Sunspots 

Danielson and Savage (1968) suggested that the missing 
flux might be carried away from the spot by horizontally propa- 
gating magnetosonic waves. Outside the spot, these become 
sound waves, and it is important to test the hypothesis cf 
whether the missing flux might be sound waves, rather than 
turbulent convection as in the preceding sub-section. 

The important quantity is the dissipation length for 
sound waves travelling through the undisturbed solar gas out- 
side a spot. Formulae for absorption of sound in a medium of 
very high thermal conductivity (such as the sun, where radiation 
ensures high conductivity even if turbulence were totally absent) 
are given by Landau and Lifshitz (1959). The result sim- 
plifies in the limit of wave frequency much greater or 

o 

much less than the characteristic frequency co c * C g /«. With 
a ■ 0.1C g ax, we have oo c * 10C g /ax. At depths of 1-10 thousand 
km below the surface, C s * 12-30 km/sec, and so with a “ 10^ cm, 
u> * (0.012-0. 03) /x. Danielson and Savage (1968) note that 
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wave periods at great depth in the sunspot are expected to be 
about 1/2 hour, i.e. u> * 0.0035. Hence with x <s#l (see previous 
sub-section), <o « cjd c , and the absorption length is 

2cj 20cj 

a ^ S CO^Q'('y-l) 0^ax(7-l) 

where 7 is the ratio of specific heats. Setting y < 5/3, x » 1, 
a * 10^ cm, C g * 20 km/sec, co * 0.0035/sec, we find 10~® cm. 

It therefore appears that low frequency sound waves could pr .pa- 
gate away from the spot ever several spot radii before dissipa- 
tion. Hence sound wave propagation away from spots is a me chan" 
ism which is worth considering as a possible explanation of th a 
missing flux in spots, 

e. Chromospheric Heating Above Spots by Alfven Waves 

The problem of chromospheric heating above sunspots in- 
volves not the horizontally propagating waves, as in the 
previous subsection, but vertically propagating waves. Since 
the field lines in the deeper layers of a spot are close to 
vertical, it is natural to consider the waves in question as 
Alfven waves. Such waves have the advantage of very long dis* 
sipation lengths in spots (Osterbrock, 1961), so that once 
they are produced, these waves can essentially propagate 
undissipated all the way to the surface layers of the spot. 

Thus the flux of Alfven waves, once they are generated, she old 
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remain constant with height. This can be proven analytically, 
using the WKB approximation (Alfven and Falthammar, i963) 
which shows that in a medium where the density and magnetic 
field varies only very slowly ever on- wavelength of an Alfven 
wave, the amplitude of the wave, b, is independent cf the field 
strength and varies as (Cowling, 1953), Teen the energy 

flux in the waves, F^-wb^V^, where is the Alfven speed, is 
independent cf density, and F, vB } the mean field strength. 

In a spot model, such as that used fer the sunspot models 
(Mullan, 1974} wi- ere B is taken to be independent of height, 
is constant . The conditions fer the WKB approximation to 
be valid are satisfied to a certain extent throughout the 
depth of the spot, but are especially well satisfied near the 
top of the spot where nrst of the Alfven wave flux is generated 
(see Figure 7 cf Mullan, 1974;. Thus the wavelength of the 
waves, (see Figure 3 of Mullan, 1974) is essentially one 
pressure scale height, H^, and so the condition for the WKB 
solution to be valid is that Hp << Hp, where Hp is the density 
scale height. It is easy ho show that 

Hp - 1 . (1-Q) V 

H P 

where *7 * d inl/d inP is a dimensionless temperature gradient, 
and Q » d £n *i/d InT is a dimensionless gradient of molecular 
weight. Due to ionization effects, Q < 0, so an upper limi. 
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on Hp/Hp is 1- < 7 . In deep layers where ionization is almost 
complete V* 0.4, and Bp exceeds H p by a factor of almost 2. 

In layers near the surface, where ionization is in progress 
| Q { increases and V approrches unity (or even exceeds it in 
seme models). Then H p becomes touch greater than H p , and WKB 
conditions are well satisfied. 

So far, our models have dealt only with condi t tens below 
the Wilson depression, i.e. at depths > 1000 km, where the field 
can reasonably be expected to be uniform. Now we must consider 
the propagation of the Alfven waves emerging froir che "surface 
of the spot" (i.e. the Wilson depression) as these waves enter 
atmospheric layers where (i) the field is no longer uniform, 
and (ii) the density varies so rapidly with height that the 
WKB approximation no longer holds. 

As regards (ii), Fikelner and Livshits (1965) claim that 

/ | /O 

the flux of Alfven waves should vary as ' . They derive 

this result from formulae which are applicable to reflection 
and transmission of Alfven waves from a discontinuity in den- 
sity (cf, e.g. Alfven and Falthammar, 1963). Thus Pikelner 
and Livshits argue that the variation of density in the visible 
layers of a spot is so rapid (H p < wavelength) that it is per- 
missible to consider the gas density as varying discontinuous ly 
from the photosphere into the chromosphere-corona. Formulae 
have been derived for this case if the field remains perperi-^cular 
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to the boundary and of constant strength across the discontinuity. 
Then the flux transmitted is 

^Trans PlP2 

^incident ( /p^ + JP2' 

where p^, p 2 are the densities on either side of the discon- 
tinuity. In Figure 6, the ratio T is plotted as a function of 
Pl/p 2 » It is seen to decrease slowly at first (p^/p 2 *tel), in 
accordance with the WKB solution that if the density varies 
infinitely slowly, the flux is constant, and all of the inci- 
dent flux must be transmitted, so that T = 1. The larger 0^/p o 

becomes, the farther the curve departs from the WKB solution 

l/2 

(I * 1). However, it is clear that to write F^~p ' , as 
Pi.kelner and Livshits have done, is too drastic, and results 
in transmission of coo little flux to the less dense medium. 

In the limit p]/p 2 >:> I* Pikelner and Livshits underestimate 
the transmitted flux by a factor of 4; and in the opposite 
linit, p^/p 2 *l> their solution does not approach smoothly to 
the WKB limit. 

As a specific example, we may consider the case of an 
Alfven wave generated at the bottom of the convection zone, 
where « 3.1 x 10 g/cm^. By the time this wave reaches the 
"top” of the spot (i.e. the Wilson depression, see Mullan, 1974/ 
where P 2 * 4.6 x 10 g/cm , it has traversed gas where p^/p2=6,7. 
In Figure 6 this point is indicated on all curves. The assumption 
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in our model of sunspots is thst T * 1, Which is probably a valid 
solution since WKB conditions are well satisfied. In the opposite 
extreme, if the density variation from bottom to top of the spot 
were encountered by the Alfven wave as a discontinuity, then 
2C% of the waves would be reflected, and X * 0.8. At most, then, 
net more than 20% of the upcoming flux should be lost, whereas 
according to Pikelner and Livshits, more than 60% of the upcom- 
ing flux should be reflected. Thus their estimate of reflected 
flux is three times larger than that predicted even in the 
extreme case of replacing the density gradient by a discontinuity. 
Since in fact most Alfven waves in the spot models are generated 
near the "tep" of the spot, pj/p 2 Is even smaller than 6,7, and 
the I curve approaches even more closely to the WKB solution, 
so that Pikelner and livshits ' estimate of reflected flux be- 
comes even more incorrect near the top of the spot. 

A more serious criticism of Pikelner and Livshits' approach, 
however, is that they claim that they are considering energy flux , 
whereas in fact they have introduced the area of a flux tube (in 
order to avoid considering field variation'*) so that their result 
actually applies to luminosity (ergs/sec) rather than flux (ergs/ 
cm /sec). When their formulae are rewritten in terms of flux, 
it is found that the flux must include the field strength, which 
certainly varies in the upper layers of a spot. However, the 
formula for T above was derived in the case of constant B, * .J 
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is therefore not applicable to the present case. It is there- 
fore not permissible to consider chromospheric heating by 
Alfven waves, and by fast and slow magnetcscnic waves, accord- 
ing to the method described by Pikeiner and Livshits, and later 
developed numerically by Marik (1967). 

It is certainly to be expected that it is permissible to 
replace a continuously variable medium by a discontinuity in 
order to study wave reflection, if the wavelength is short com- 
pared with the density scale height. Justification for this 
expectation is provided, at least in the gas -dynamic case, by 
the work of Chisnell (1955) who shewed that when a shock impinges 
cn a region where the density varies smoothly from pj to p 2 , the 
strength cf both reflected and transmitted waves agrees well 
with the results obtained if the shock impinges on a discontin- 
uity where the density rises abruptly from tc p 2 ^ A shock 
is however a more extreme case than we are considering here, 
and the principal problem associated with Alfven wave propaga- 

t 

t 

ticn in a spot is this: the wavelength of the wave, X; is 
initially (at great depths) small compared with the local 
density scale height, X « H^, but in propagating upwards the 
Alfven speed grows, and the density scale height falls, until 
eventually near the photosphere (ZssO, on the scale where the 
Wilson depression is Z^-1000 km) X becomes larger than Hp, 

In the topmost layers cf the upper photosphere and lew chroma 


1 


i 



I 


1 

l 


i 

\ 


1 


f 




{ 


J 


-30- 


sphere, assuming that tke wave frequency remains equal to that 
at excitation (period **10^ - 10^ seconds), X eventually becomes 
extremely large compared with Hp» It appears that only a num- 
erical solution can cope with conditions where X/Hp assumes 
values which are both less than unity and greater than unity. 

The theoretical picture becomes more complicated if the 

field is not perpendicular to the discontinuity. All three 

hydromagnetic wave modes are coupled together, especially if 

the incident waves move nearly parallel to the field lines 

(Stein, 1971). The solid angle centered on the field lines 

within which coupling of modes is strong is proportional to 
2 

(Ap/p) where Ap is the discontinuity in density at the boundary. 
In the limit of continuous variation of density, the coupling 
goes to zero unless the direction of the magnetic field ro- 
tates with altitude. In the latter case Alfven waves propagat- 
ing through a stratified medium couple with the fast -mode waves 
(Frisch, 1964). 

The propagation of hydromagnetic waves in a stratified 
medium in absence of mode-coupling (i.e. in the absence of field 
rotation) was discussed by Ferraro (1954), and later 
by Ferraro and Plumpton (1958). They assumed 
a uniform field, and applied the boundary condition that the 
wave amplitude must go to zero at large altitudes. The solution 
depends on the ratio of X to Hp. When X « Hp, the amplit~n= 
of velocity in the wave varies as predicted by the WKB approx- 
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imation v « at high densities. Eventually, however, the 

decreasing density causes X to exceeo Hp, and then the velocity 

1/2 

amplitude decays rapidly to zero, v « p ’ . Thus even in a 

uniform field, the flux of energy would decay proportional to 
3/2 

p in the upper atmosphere. The wave distorts as it propa- 

* 

gates upwards, in the sense that gravity inhibits vertical 
motions and causes the wave to become polarized horizontally. 

Hollweg (1972) criticized the boundary condition used by 
Ferraro and Plump ton, and proposed that a more correct boundary 
condition is to require that there are no waves coming in from 
infinity. Hcllweg retained the case of constant vertical field, 
but extended their analysis to a stratified medium with twc 
scale heights, one small at low altitudes, the other large at 
great altitudes. This is a good approximation to the solar 
chromosphere. The solution found by Kollweg behaves as expected 
in the WKB limit, namely v « However, in the non-WKB 

limit, i/Hp >> 1, the amplitude continues to increase with 
height, although mere slowly than in the WKB case: v* log (x/Hp). 
The normal hydromagnetic relation between velocity amplitude v 
and magnetic amplitude b does not hold in the non-WKB case, 
and b becomes constant in the upper atmosphere. 

It is a feature of ncn-WKB propagation that reflection 
can occur even if the density is continuous. A discontinuity 
in scale height serves to reflect a non-WKB wave. The flux 
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1 . 


leaking upwards into the chromosphere depends on how much is 
reflected at the discontinuity in scale height. For the special 
case considered by Hollweg, only about 1/6 of the upward moving 
flux can leak through to the chromosphere. We have extended 
Hollweg' s results to the general case of the discontinuity in 
scale height occurring at an arbitrary altitude. The transmis- 
sion coefficient T can be obtained in the limit of wavelength 
much larger than scale height by using the asymptotic expansions 
cf the Jq and Yq Bessel functions: 

T » , My 

(0+ yy" + JL. (y/ny-0in0) z 
x 2 

where 0 * 2 Zq^o)/Vq 2 and y * ^ z 02^ z 01* Here, 1 * 03 *- s t ^ ie wave 
frequency, Z ^ is the scale height near the photosphere, z 02^ >z qi^ 
is the scale height in the chromosphere, and Vq 2 is the Alfven 
speed at the altitude where the scale height is discontinuous. 

If the discontinuity in scale height is shifted up in the chromo- 
sphere to a density 100 times smaller than that considered by 
Hollweg, we find T = 0.05, which is smaller than the value 
obtained by Hollweg by a factor of only 3 or 4. Thus the 
leakage of non-WKB waves into the chromosphere is not sensitive 
to the density at which the chromospheric temperature rise 
sets in. The altitude of the temperature minimum can alter by 
several hundred kilometers without essentially changing the 
transmission coefficient, at least in the case of constant 
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magnetic field. 

In an attempt to extend Hollweg’s analysis to the case of 
a non-uniform field,, we consider Hollweg's equation for the wave 
amplitude v as a function ot the modified height variable ?*e” Z ^ 2Z 0: 

I 12 

• v" + f* 1 v* + a* v m 0 (K-2) 

i 

! 

j where prime denotes differentation with respect to ?, and 

a 1 - 4Z o“ 2 /v o- In Hollweg's case, Vq is the Alfven speed at 

height Z * 0, and is therefore proportional to Bq, the assumed 

constant field. The case of varying field can be represented 
2 

by allowing a to vary with altitude. Since B is expected 

to decrease at great heights, « 2 must increase as Z increases 
2 

i.e. Q must increase as IT decreases. Suppose we write then 
«2 * a 2 + 

where is the value of a deep in the photosphere. Then 
equation (H-2) becomes 

v" + v' y _1 + (I+jj2/y 2 ) V = 0 

where y « and primes denote differentiation with respect 
to y. The above equation has a solution which can be considered 
as a Bessel function of imaginary order. The solution is given 
by Boole (1844): 

v * cos (pjny) [ 1 + y 2 + a^ y^ + . . » J 
where coefficients a^ etc. reduce to those appropriate for 
the zeroth order Bessel function Jq(v) in the limit p, * 0. There 
is a corresponding series solution, singular at y ■ 0, corresponding 
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to Yq(y). The numerical values of «q can be obtained approxi- 
mately by setting 5 * 1, Zq * 185 km, a> - 0.02 sec**, Vq ■ 

10® cm sec * (appropriate for a sunspot with B » 3000 gauss, 
p » 7 x 10 gm cm J ) . Then «q * 0.74. At the top of the 

atmosphere, say Z « 2000 km, where B has decreased to 300 gauss, 

2 2 
a must have increased by a factor of 100, so that u » 

2 , . -4 

99 exp (-2000/185), or p.»10 . Then the modification 
of the Jq solution for v due to non-zero p. is the factor 
cos 10^ Jin «q -(Z/Zq) . As long as we are interested 

only in wave propagation within a few tens of scale heights 
distance fiom Z * 0, this factor cannot depart from unity by 
more than a few parts in 10^. This suggests that as long as 
the field varies with altitude in a way which is much slower 
than the variation of density, then the full solution of the 
wave equation should approximate the (Jq, Yq) solution obtained 
by Hollweg. However, the best test for this is to obtain a 
numerical solution. 

We have attempted to derive a numerical solution as 
follows. The usual hydromagnetic equations are written in 
the form (Deift and Goertz, 1973): 
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¥ * <P*0 - 0 

where f is the potential including (if necessary) gravitational 
and centrifugal terms. The field of the spot is written in 
cylindrical cooruinat* s as Bq * (B r , 0 S B z ). The field of the 
Alfven wave, and its associated velocity amplitude, are written 
as (b r , 0, b^) and (^ r , 0, v z ) respectively. In the absence of 
waves, • Bq « 0, and so to satisfy 7 • B * 0 in the presence 
of the waves, we have V • b * 0. It is assumed that in the 
absence of waves, the field satisfies 

7(SL+,tf). So -VBq 
- 8 » 

This is the assumption of magnetohydrostatic equilibrium, includ- 
ing curvature forces, and it presupposes that the spot is a 
stable entity. It is assumed that in the perturbed state we 
may keep terms of first order only. Then since variations of 

pressure in an Alfven wave are of second order (Alfven and 

" / 

Falthammar, 1963) we neglect 5p. Also, in an Alfven wave 

o 

Bq • b * 0, and so to first order B remains unchanged. Assum- 
ing that variations in 0' are also of second order, the perturbed 
version of equation (2) is 

1 

^ ) first order terms 

combining this with the perturbed version of equation (1), and 
noting that b r /b z ■ ~B Z /B r in a wave where bg is assumed to remain 
identically zero, we obtain the following equation for b r : 




where the wave has been assumed to be sinusoidal with frequency 
cd, and the quantities U and V are the following: 


U - B 4 /B 2 

Z 

- iL. (, ^ £ r ^ 3 B r _ ^r ^ B z _ B r ^ B z i 

r B z B z Jr B z 2 z B z Jr Jz ( 

and B 2 « B 2 + B 2 . 

r z 

As background fields, we have taken that suggested by Yun (1970): 
B z (r,z) - D(o) ? 2 
B r (r $ z,) * 

0 ( 0 ") *» (<I>/ ?r ) exp ( -a 2 ) 
a * r?(z) 

Uz> - ( ir B z (0,z)/C , ) 1/2 

$ » magnetic flux through the spot # 

Once the vertical field along the axis is specified, B z (0,z), 
these equations define the magnetic field at all r and z. 

Along the axis we have taken B z (0 : ,0) * 3000 gauss, and a grad- 
ient of *1 gauss/km so that the spot field is no longer vertical 
at altitudes of 3000 kc above the center of the spot. 
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To solve for the amplitude b r , two boundary conditions are 
necessary. At the moment, the two we have chosen are the fol- 

O 

lowing: (i) b r (0,0) *• 10 gauss; this ensures that the local 
flux of Alfven waves at the center of the spot at z » 0 is 
approximately equal to that generated in our sunspot model with 
B « 3000 gauss (Mullar., 1974); the numerical value of b r , 
though, is simply a scaling factor as far as the solution is 
concerned; (ii) assume the umbra is uniformly dark, then the 
flux of Alfven waves must be constant at all radii at z = 0 
(in this integration, z * 0 is located at the level of the. 

Wilson depression, at which we presume we know the flux of 
Alfven waves entering from below); thus 

b 2 B * constant =■ b 2 (0, 0)B Z (0,0) = 3 x 10^ 

This leads to 

b r (0 3 0) exp (o ,2 /2) 

b r (r,0) « 

Hence «?b r /Pr is known at z ** 0, and '<7 • b » 0 then leads to an 
expression for db^l&z at all radii in the umbra at z s 0, At 
each radial position in the umbra, then, equation (3) can be 
integrated upwards numerically, once the variation of density 
with height has been specified. We assume (see Mullan, 1974) 
that the density in the spot is equal to the density in the 
undisturbed sun, so that the density could ideally be read off 
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from, say, the BCA model (Glngerlch and de Jager, 1968), How- 
ever, numerical differentiation of this model leads to a non- 
monotonic variation of density scale height (cf. Jordan, 1970), 
and since it is preferable to keep all functions varying as 
smoothly as possible, we have decided to choose Hp initially 
to have a smooth variation with height, following as closely 
as possible the BCA variation for those regions of the atm^s- 

\ 

phere where BCA is relevant. In the deeper layers, near Z «■ 0, 
we have used interpolated H p values between the deepest BCA 
I value and the value in our sunspot model at the Wilson depres- 

! sion. Then choosing the density to be equal to BCA density at 

the level of minimum H n (p = 1.^3 x 10”^ gm/cm^ at H = 107 km), 

l p p 

! we derive the model as in the table; 
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TASLE 2 
pit m/ca?) 


0 

8.79 x 10“ 7 

1900 

150 

7.93 x 10’ 7 

1190 

300 

6.83 x 10° 7 

832 

450 

5.53 x 1C* 7 

584 

000 

4.09 x 1CT 7 

409 

750 

2.65 x 10“ 7 

286 

900 

1.43 x 10* 7 

200 

1050 

5.52 x 10°® 

115 

1200 

1.43 x 10“ § 

107 

1350 

3.54 x 10“ 9 

108 

1500 

8.89 x 1C‘ 10 

109 

1650 

2.27 x 10“ 10 

111 

1800 

5.96 x 10“ U 

113 

1950 

1.61 x 10“ U 

116 

2100 

4.51 x 10“ 12 

120 

2250 

1.33 x 10* 12 

125 

2400 

4.12 x 10" 13 

132 

2550 

1.40 x 10-1 3 

145 

2700 

5.39 x 10* 14 

170 

2850 

2.39 x 10* 14 

200 
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Tfce height scale is such that z » 796 corresponds tc "Height * 0" 
in the BCA model. The density in the table follows the density 
in the BCA model in the chromosphere very well, departing from 
it by no more than 25% even at the topmost heights. The deeper 
layers do not agree as well with the run of density derived in 
our sunspot models. Thus at z *• 0, the density in Table 2 
exceeds that in our sunspot model at the Wilson depression 
by 85%. However, this is net excessive, and the solution 
using Table 2 should preserve the major features of the solu- 
tion using a more detailed variation of density with height. 

The period chosen for the waves is 300 seconds. This 
is the period cf waves excited at the Wilson depression, in 
cur spot model if we set the teried equal tc the local cell 
height (Hp) divided by the convective velocity (v ). It takes 
gas a period of time cf this order tc circulate around a cell 
and thereby generate an Aifvett wave according to the scheme of 
Figure 3 of Mullan (1974). In a spot with B * 3000 gauss, we 
find Hp s ^220 km, v c &0.75 km/sec. ^therefore, #p/ v c “ 294 
seconds. In spots with B = 2000 and 4000 gauss, we find Hp/v * 
240 and 35C seconds respectively. A period of 300 sec is long 
enough thac internal gravity wave excitation in the atmosphere 
is expected to reflect some of the Alfven waves back down into 
the spot. Most of the waves should be transmitted, however 
since I * 70% (Danielson and Savage, 1968, especially Figur.- . 
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in their paper). We therefore exp® 1 ' the wave flux upwards to 
be reduced by about 30%. At first , the numeric*.! solution for 
the flux on the axis cf the spot dees indeed show a decrease 
with height, dropping off tc 65% cf incident at a height of 
600 km above z * 0, i.e. at a dept ', of about 240 km below the 
undisturbed photosphere. Above z * 60C km, however, the flux 
begins to increase with height (see Figure 7;. 

Away from the axis of the spot, where the field direction j 

i 

is not vertical, the vertical flux is smaller than on axis (at [ 

i 

z*0). The flux at these positions, however, alsc shows an 
increase cf flux with height (see Figure 7). 

These flux increases are net physical y acceptable, buc 

£ 

it is not yet clear what is causing the increases. Numerical 
errors may be involved, cr perhaps seccnd order terms need to 
be retained in the perturbation equations. Another possibility 
is that we have encountered a resonance peak of the type discussed 
by Hollweg (1972, especially his Figure 2). The principal peak 
occurs when * 2 Zq^/Vq^ equals the first zero cf the Bessel 
function Jq. In our case, with cc * 0.02, Vqi * 10 km/sec, we 
find Jq(g) * 0 if Zq^ 500 km. At such a peak, the flux may 
rise tc ten times (or more) the WKB flux. Another possibility 
is that mode-mode coupling becomes important when the wave- 
length of the wave (* * 220 km * pressure scale height at the 
Wilson depression) approaches, and becomes less than, H p . 
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According to Table 2 above, this should occur at z * 800 km. 

For any reasonable wave period, the Alfven wavelength A 
inevitably exceeds the density scale height K p in the upper 
atmosphere, although the wavelength is less than the density 
scale height deep below the surfs te at (say) the level of the 
Wilson depression. In an attempt to make the transition from 
^Hp to smoother, we have computed model atmospheres for 

umbrae using the computational method of Mihalas (1967). We 
have assumed that temperatures in the umbra differ from these 
in the Harvard -Smithsonian Reference Atmosphere (HSRA) by 
A? * constant (where 6 * 5040/ T). (This is equivalent to 
assuming radiative equilibrium in the umbra.) We have com- 
puted models with A 0 * 0.1 (0.1) 0.6 using the opacity tables 
of Bode (1965). Some of these me dels are shown in Figure 8, 
and for comparison, the model of Kneer (1972) is alsc shown, 
Kneer'r model is indistinguishable from the model A0 * 0.4 
high in the atmosphere, where radiative equilibrium must 
exist. Deeper in the umbra, the departure from radiative 
equilibrium is due to the onset of convection. 

For purposes cf our Alfven wave computations, the 
den' le heights are of interest (see Figure 9). In the 

solar moael s the presence of a density inversion causes H 

r 

to be a very irregular function of optical depth. Such behav- 
ior would greatly complicate the investigation cf alfven wav- 
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propagaticn. Fortunately, in the uzabral models having A© > 0.3, 
the irregularities disappear, and the density scale height now 
varies smoothly over the entire range of optical depths consid- 
ered here. This is not to say however that irregularities 
would not be encountered at greater depths. In the solar model, 
the inversion coincides with the onset of hydrogen ionization, 
and in the umbral models, this is postponed to greater depths. 

Unfortunately, no convergent solution has yet been found 
for the Alfven wave flux. The above equation (3) includes 
only terms of first order in the wave amplitude. Second order 
terms have now been added, but the program is new no much more 
complicated that debugging has not yet been completed -satisfac- 
torily. At present, the solutions have not yet been obtained. 

Rotation of the plane of polarization of the wave must also be 
considered. 

f . Magnetic Convection in the Sun and Solar Neutrinos 

It was suggested by Schatten (19 '3) that magnetic fields 
inside the Sun might rise to the surface by buoyancy from the 
deep interior, thereby contributing to the transport of sclar 
energy flux. He speculated that this might have a sufficiently 
large effect on the temperature structure of the sun that the 
neutrino flux emitted by the sun would become small enough to 
agree with observations. We have shown that this is very un- 
likely, for the velocities of upward motion of mag'.ietic flux 
ropes in the radiative interior can be shown tc be cf order 
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I 0.01 cm/sec or less. Such small velocities cannot change the 

temperature structure by mere than one part in 10^. In particu- 
lar, the central temperature will not change by mere than a 
small fraction of 1%. Thus the effect of magnetic fields 
inside the sun can in all probability not contribute signifi- 
cantly to resolving the neutrino controversy. 

g. Models of Starspots cn Flare Stars 

In last year's annual report, we described a model devel- 
oped for sunspots based on the assumption that the missing flux 
was in the form of Alfven waves (Mullan, 1974). The waves were 
assumed to be generated by convective gas motions, which occur 
with reduced efficiency in the presence cf a vertical magnetic 
field. A quantitative estimate of the reduction in convective 
efficiency is incorporated by means of a physically meaningful 
modification of the cellular model of turbulent heat convection 

it 

due to Opik (1950; also Mullan, 1971). The model had led to 
results for sunspots which were in agreement with a variety of 
observations . 

This model has now been extended to the case of cool spots 
on stars of late spectral type on the lower main sequence, and 
the results agree with observations very well. Further applica- 
tion cf the starspot model to another class of stars will be 
described in the following section. Here we confine our r err. arks 
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to the flare stars. 

The model in the case of sunspots had indicated that for 
field strengths less than a critical value B c , no spot model was 
possible. This lower limit arises because for fields less than 
B c> the field is strong enough tc reduce the efficiency cf con- 
vection appreciably, and yet not strong enough tc compensate j 

i 

•ior the reduced convective flux by emitting sufficient Alfven j 

waves. Thus, flux constancy cannot be achieved unless B B c ^ 

1200 gauss in the sun. In a red dwarf of spectral type ' 

MO vl e ff * 4000 C K> , B c was found to be much larger, of order 

i, 

B c * 10-15 kilcgauss. Hence starspots on rod dwarfs require 

much stronger fields than sunspots. However, the fields cannot 

become too strong, for then the Wilson depression exceeds 1000 

km and the spot would be washed out rapidly by radiative heat- \ 

i 

ing across the walls of the spot. This sets an upper limit of 
30-40 kilogauss cn the field strengths in starspots. A typical 
surface field strength of a starspot cn an M0 star (with 
masses 0.6 M ) is therefore about 20 kilogauss. 

Independent estimates of field strength in such a star can 
be derived from flare energies and from the general properties 
of the variation of pressure with depth. Ihese both ie.id tc 
estimates of field strengths in agreement with the above estimate. 

More recently, in considering the evolution of stellar dynamos, 
we have in a completely independent way verified that field 
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of order 20 kilogauss are to be expected in starspots on stars 
with M - 0.6 M 0 « 

With a surface field of 2 x 10^ gauss, it is necessary to 
incorporate a field gradient so that the field becomes stronger 
inside the star, and makes the starspot as deep as the convection 
zone. It is a postulate of our model (Mullan 1973) that a spot 
is a special kind of convection cell which penetrates to the 
base of the convection zone, and has a diameter, D, which is 
simply proportional to the depth of the convective zone, K, The 
constant of proportionality, D/H, can be estimated for poly tropic 
models, and in the case of the sun, D/H varies from 2 to about 
7, with preferred values apparently at 2.9 and 5.4 (cf. Bumba, 

Ran zinger and Suda, 1973, especially their Figure 7, and note 
that we have taken the depth of the solar convection zone to 
be 10^ km, following Mullan, 1971). This leads us then to 
exract that in stars of later spectral types, where the convec- 
tion zone occupies an increasingly large fraction of the radius 
of the star (see Figure 10 ) f the fraction of the surface area 
occupied by a spot must increase. In fact, in the cases of 
the stars YY Gem, CC Sri and BY Dra, the sizes of the starspots 
(Bopp and Evans, 1973; Bopp, 1974) are consistent with our 
hypothesis if D/H = 3, 3, and 5 respectively. These are close 
to the preferred solar ratios of 2.9 and 5.4, and this gives us 
some confidence that our postulate of D proportional to R it 
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acceptable. Further confirmation of this postulate will appear 
in the following section. 

Once the field strength is known at the surface, and a 
gradient has been chosen, the computation of a spot model proceeds 
as for the sunspot case, except that the condition of equal gas 
density inside and outside the spot has been relaxed. Results 
are indicated in Figure 11. The effective temperature of the 
spot at the surface turns out to be very low, about 1600 C K. 

Bopp and Evans (1973) had suggested that the effective tempera- 
ture in spots observed by them might be about 2000°K, but this 
number is extremely uncertain, and the good agreement with the 
umbral surface effective temperature deduced here must be con- 
sidered coincidental, 

h. Starspots on the Primaries of Contact Binary Systems 

Further application of our starspot models is new possible 
in the case of W Ursae Majoris type binaries. These binaries 
are considered to be encompassed by a common convective envelope 
which has the result of equalizing the entropy in each compon- 
ent's convection zone (Lucy, 1968) with the result that effec- 
tive temperatures of both primary and secondary components 
become equal. Observations cf hotter W *JMa systems (the so- 
called A-type systems) indicate that the two components do 
indeed have closely equal effective temperatures. On the otr =*.r 
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‘ i 

' hand, the cooler systems (called W-types) are obeerved to have 

temperature differences between primary and secondary, AT, 
which are of order several hundred degrees. In terms of the 
parameter X * AT/T, the observations suggest that X ranges from 

I 0.02 to 0.11-0.13 in different W UMa systems. 

. j 

S It is possible to interpret these X-values in terms of 

! 

| starspots on the primary if the spots have diameters proportional 

j 

j to the depth of the convection zone. The preferred values D/H ® 

j 2.9 and 5.4 lead to theoretical values of X as a functi^ of 

i stellar B-V color as shown in Figure 12. Points denote observed 

values. It seems that our starspot postulate leads to a very 
plausible interpretation of the photometric peculiarities of the 
W-type systems. 

The presence of spots causes distortions of the light 
curves, and we have obtained quantitative agreement (to accuracies 
of 0?02) between predicted and observed distortions in four 
systems. 

The W UMa systems are especially suited for stellar 
dynamo activity on account of the simultaneous occurrence of 

i 

deep convection zones and rapid rotation. The latter is en- 
forced by tidal effects. 
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i. Implications of Starspots on Red Dwarfs 

If the fields in starspots on red dwarfs are indeed as 
large as 20 kilogauss, then fields in active regions around such 
spots must also be much larger than in the sun. This leads us 
to expect that hydromagnetic heating of chromospheres in magnetic 
red dwarfs must be about two orders of magnitude greater than 
in the sun. We believe that this explains why chromospheres of 
dMe stars are 10-100 times denser than in the sun. 

Such dense chromosphere-coronas also lead us to expect 
that analogous to the Moreton waves which are emitted following 
certain large solar flares , and which may trigger sympathetic 
flares in distant active regions, stellar flares also emit 
Moreton waves which travel with speeds roughly equal to those 
observed in the sun. 

Starspots have effective temperatures which are so cool 
that grains in all likelihood can condense. If dMe stars have 
dense chromospheres because of strong magnetic activity, then 
grains should form in dMe stars following starspot activity. 
However, dM stars presumably have no magnetic activity, and so 
grains would not be formed. We believe that it is grains in 
dMe circumstellar shells which causes the slight systematic 
infrared excess of dMe stars relative to dM stars. 

A point of major interest in stellar chromospheres is the 
question of possible cyclic activity. No such cycles have yet 
been definitely determined. We believe that this is due to the 
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long lifetime of starspots on these stars. Vogt (1973) has found 

that a spot on BT Dra has survived as long as eight years. If 

starspots survive longer than the period of the magnetic cycle, 

it will be very difficult to discover cyclic behavior. On 

the b^sis of a numerical model of the evolution of stellar 

dynamos in red dwirfs, we believe that the cycle time in red 

dwarfs may in fact be quite short, or order 1-3 years, and 

therefore much shorter than the lifetime of starspots. 

The release of energy in stellar flares can be on a much 

34 

grander scale than in solar flares. Thus 7 x 10 ergs of 
radiative energy was released in a flare of YZ CMi (Kunkel, 

1969). This exceeds the largest solar flare by a factor of 

3 

almost 10 . It is hoped that an understanding of the relation 
of starspots and stellar flare may help to understand the rela- 
tion between sunspots and solar flares, and vice versa. 

IV. Personnel Activities 
Paper Presentations and Meeting Attendance 

"Observations of the Sun", by D.J. Mullan, colloquium talk at 
the Physics Department, Lycoming College, Williamsport, Pennsyl- 
vania, November 23, 1973. 

"Starspots on Red Dwarfs", by D.J. Mullan, colloquium talk at 
the Astronomy Department University of Texas at Austin, April 
22, 1974. 

"Magnetic Fields in the Sun", by D.J. Mullan, invited paper at 
a symposium commemorating the 50th Anniversary of Bartel 
Research Foundation, May 11, 1974. 
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"Magnetic Fields on Spotted Red Dwarfs”, by D.J. Mullan, paper 
presented at the 143rd Meeting of the American Astronomical 
Society in Rochester, New York, August 23, 1974, 


Papers Published 

”Flare Triggering by Coherent Oscillations”, by D.J. Mullan, 
Astrophysical Journal 185 , 353, 1973. 

"Sunspots, Supergranules, and the Depth of the Solar Convection 
Zone”, by D.J. Mullan, Astrophysical Journal 186, 1059, 1973. 

"Sunspot Models with Alfven Wave Emission”, by D.J. Mullan, 
Astrophysical Journal, 187 , 621, 1974. 

"Starspots on Flare Stars”, by D.J. Mullan, Astrophysical 
Journal, 192, 149, 1974. 

"Correction of Sunspot Intensities for Scattered Light”, by 
D.J. Mullan, Solar Physics 32, 65, 1973. 

"Conments on Papers by P.R. Wilson Concerning Sunspots”, by 
D.Jo Mullan, Solar Physics 32, 441, 1973. 

"Fast Rotation of Metal Poor Stars", by D.J. Mullan, Astronomy 
and Astrophysics 2J _ , 379, 1973. 

'Magnetic Fields on Spotted Red Dwarfs", by D.J. Mullan, Bulletin 
of the American Astronomical Society, 6, 333, 1974. 


Papers in Press 

"Magnetic Fields in the Sun”, by D.J. Mullan, Journal of the 
Franklin Institute. 

"Magnetic Convection in the Sun", by D.J. Mullan, Solar Physics. 
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Papers Submitted for Publication 

"Magnetic Fields and Dense Chromospheres in dMe Stars", by 
D.J. Mu 1 lan, Astrophysical Journal, 

"On the Possibility of Magnetic Starspots on the Primary Compon- 
ents of W Ursae Majoris Type Binaries", by D.J. Mu 1 lan, Astro- 
physical journal. 

Personnel Changes 

Dr. Arne A. Wyller left Bartol to take up responsibilities 
as Research Professor and Director of the Swedish / strophysical 
Observatory at Anacapri, Italy. 

During July and August 1974, an undergraduate summer 
student, John Oliensis, of Yale University, contributed to the 
tneoretical work on Alfven wave propagation. 

Dr, D.J. Mullan has continued full-time work on both the 
observations and theoretical aspects of sunspots and starspots. 
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Captlcns for Figures 

Fig. 1 Idealized spectrum of spot near the sodium D lines, 

indicating the intensities involved in defining a two- 
dimensional classific. ticn of sunspot atmospheres. 0 
denotes continuum windows in sunspot spectra, following 

t« 

Wohl (1970). I A is the interpolated continuum inten- 
sity in the spot at 5893. 4A. lg is the actual measured 
intensity in the spot at 5893.4 a. The classification 
suggested here is r g p 0t * Ig (spot)/l A (spot) and 
Rspot * I A (spot)/l A (sun). 

Fig. 2 Two-dimensional classification of the line-snd continuum- 
forming regions in sunspot atmospheres. See legend of 
Figure 1 for definitions of r g p 0t and Rsp 0t; . Observations 
are characterized by lines, indicating the date of obser- 
vation, the umbral area in millionths of a solar hemis- 
phere, and the minimum Z correction factor for scattered 
light. Note that errors due to scattered light lead to 
sloping lines in this figure, rather than error boxes. 
Theoretical model predictions at disk center are indicated 
by enclosed letters (K * Kneer (1972); SW * Stellmacher 
and Wiehr (1972); Y » Yun (1971); YI, YII » revised models 
of Yun). Two observations with large r g p 0t; ( > 0,9; are 
near the limb, where theoretical models predict iar^-r 
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Fig. 3 Broad-band intensity of spots expressed in terms of the 
disk intensity. Four filters were used on October 24, 
1973, seven on October 26, 1972. Solid lines are 
theoretical model predictions cf continuum intensity. 

The large uncertainty at 4200A is due to uncertainties 
in correction for line blanketing in the spot. It is 
unlikely that errors in observations and reduction can 
account for the large differences between the two spots , 
The larger umbra appears to be definitely darker than 
the smaller umbra. 


Fig. 4 Dots represent counting rates per second U as a function 
cf position as the telescope is stepped in 2" steps north 
of the center of the umbra. At the center cf the umbra 
(point labelled C) the counting rate is 206 thousand 
photons /sec. The line ED is an eye -fitted line through 
the observations, indicating the decrease in intensity 
away from the spot. The undisturbed disk intensity 
corresponds to a counting rate of 470 thousand photons/ 
sec, although limb darkening needs to be taken into 
account across the scan. The missing flux in the spot 
may have reappeared as the excess (denoted by line ED) 
of intensity above the undisturbed intensity (denoted 
by BD'; . 
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Fig. 5 Function G(z,p) appearing in the solution for heat con- 
duction outside a cylinder of radius a. Curves have been 
drawn through values tabulated by Ingersoll et al_. <, 1950) 
for the case when a constant heat flux emerges across 
the curved surface of the cylinder. At a given time, 
excess temperature at radius r outside the cylinder is 
proportional to G evaluated at p * r/a along a vertical 
line of constant z. At a given position outside the 
cylinder (constant p), the excess temperature grows 
asymptotically as leg (time). 

Fig. 6 T * Ratio of transmitted flux to incident flux when 

Alfven waves impinge on a discontinuity where density 

decreases from tc p 2 « Solid curve is the function 

»• 

derived by (e.g„. Alfven and Faithammar i.1965; Ibe 
dash -dot curve is the approximation used by rikeiner 
and Livshits <,1965) for Alfven waves; the dotted line 
is their approximation fer fast-mode waves. In the WKB 
approximation (wavelength much less than density scale 
height) the flux of waves transmitted through a regicn 
where density declines gradually from to p 2 remains 
equal to unity (dashed curve) . The arrows denote 
special case of waves propagating between the bottom 
and the top of the solar convection zone, 
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Fig. 7 Flux of Alfven waves as function of height above Wilson 
depression. Solution is obtained by numerical integra- 
tion of equation (3) at fixed radius, assuming that the 
field geometry in the spot is that of Yun (1970). 

Numerical errors, or incorrect approximations , or both, 
may be responsible for the unphysical behavior of the 
solution in those regions where the Alfven wave flux 
exceeds the total solar flux. 

Fig. 8 Temperature structure of umbral models computed accord- 
ing to the method of Mihalas (1967). The temperatures 
in the umbra differ from the solar temperatures at 
equal optical depth by an amount A 9 * constant. The 
umbral model atmosphere of Kneer (1972) is shewn fer 
comparison. 

Fig. 9 Density scale height as a function of depth in the 

umbral atmosphere models. Irregularities in the solar 
model are caused by a density inversion in the HSRA model. 

Fig. 10 Depth of convection zones in stars on the lower main 
sequence. Depth is expressed in terms of the stellar 
radius; abscissa is mass in units of the solar mass. 
Results for M/M 0 ^ 0.4 are taken from Copeland et al_> 
(1970). Completely convective stars set in at M'M *0.3 
(Copeland et al. ) or at M/M 0 = 0.20 (Straka 19?' 1 », 
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Fig. 11 Physical variables in a model of a starspot on YY Gem 

as a function of depth. The depth is expressed both in 

terms cf the gas pressure in the undisturbed stellar 

model and in terms of the depth Z(cm). B is magnetic 

field strength; D/H is the local ratio of convective 

cell diameter to cell depth; T g is the local effective 

temperature in the spot, defined such that the local sum 

of radiative plus convective fluxes equals o t\ where 

£ 

a is the Stefan-Boltzmann constant. Seme of the fluctua- 
tions in T g are due to numerical errors when the tempera- 
ture gradient is very close to adiabatic. NE denotes 
the range cf depths in the model where numerical errors 
are so severe that there occurs an apparent suppression 
of convection altogether. This is net considered to 
be physically meaningful. However, the decrease in T e 
near the surface is real. 

Fig. 12 Effect of a starspot on effective temperature of main 
sequence stars. The curves are theoretical estimates, 

(jl, of the decrease in the mean effective temperature 
of the visible disk of the primary component when a 
starspot is fully visible on the disk. Spot diameter 
D is assumed to be proportional to depth of convection 
zone H. Numerical values of H (Copeland et al, 1970) 
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are obtained here by assuming that B-V is convertible 
to T fi using Johnson's (1966) tables. Observed values 
of X * AT/T for individual W UMa systems are taken 
from Table 2 of Rucinski (1974). Colors are taken 
from Mauder (•; 1972) and Koch (o; 1974). Temperature 
of spot varies with stellar temperature (see text), 
lumbers beside tick marks on each curve are magnitudes 
of fading in visual light caused by the spct if it 
existed on a single star. In a close binary, the 
total light will fade by a smaller amount. 
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